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CHAPTER 1 


INTRODUCTION 

The accretion of ice on the lifting or propulsive 
surfaces of an aircraft is a phenomenon which can have severe, 
if not disastrous, consequences if it is not dealt with 
properly and in a timely fashion. In many cases, this 
accretion may produce a serious reduction in lift and Increase 
in drag, requiring more power than is available and resulting 
in a serious departure from the mission profile. Although 
> measures and devices exist which generally prevent this worst- 
case scenario from happening, ice accretion on airfoils is 
still a common occurrence in certain atmospheric conditions 
and as such merits further study because of the significant 
negative effect which it has on airfoil and aircraft 
performance. It is the objective of this thesis to provide a 
single computer code which will accurately predict the degree 
of performance degradation, primarily with respect to drag 
-' increment, which will result from the exposure of an airfoil 
to some set of atmospheric and flight conditions conducive to 
ice formation. 

Two basic types of ice are found to occur as a result of 
:• exposing an airfoil in forward motion to supercooled water 
droplets in a sub-freezing environment (Fig. 1), The type of 
ice which will form may be determined by a variety of factors 
including freestream velocity, liquid water content of the 
cloud, droplet size distribution, and freestream temperature. 
Based upon an investigation of the data currently available. 
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the two latter factors appear to be the most Influential in 
determining the resulting ice type. The first kind of ice, 
known as rime ice, occurs at relatively low velocities, low 
liquid water content values (typically 0.5 to 1.0 grama per 
cubic meter), and temperatures well below freezing. Due 
primarily to the very cold temperatures associated with rime 
ice formation, the droplets tend to freeze on impact to form a 
fairly smooth addition to the leading edge of the airfoil. At 
some point with respect to the combination of flight and 
‘atmospheric conditions present, ice type transitions from’ rime 
to glaze ice. Not only is the crystalline structure of these 
ice types different, but glaze ice is also produced by higher 
freestream velocities, higher" liquid water content values (on 
the order of 1.5 to 3.0 grams per cubic meter), and 
temperatures 1 near, but below, freezing. Again, due chiefly to 
"the warmer temperatures involved, the water droplets impacting 
the surface tend not to freeze on impact but rather to strike 
the airfoil or existing ice formation and run back somewhat in 
a chordwiae direction before freezing. 

Both types of ice cause significant performance 
degradation of J the airfoil. When compared to the aerodynamics 
of an airfoil in the clean, or non-iced, ; configuration, the 
iced airfoil will virtually always exhibit a decrease in lift, 
an increase in drag, and a change in pitching moment which 
will yield a detrimental effect on aircraft stability. 
Examination of the available icing data base indicates that 
the degree of performance degradation, specifically with 
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respect to drag increment, is generally more severe for the 
glaze ice condition. This may be attributed to the more 
drastic alteration of the leading edge of the airfoil as a 
result of glaze ice formation, as compared to the typically 
smoother and less obtrusive rime ice formation. Increases in 
drag coefficient in excess of 100* or more are not unusual for 
either type of accretion. 

Both rime and glaze ice act to negatively influence 
airfoil performance by physically altering the shape of the 
airfoil, thereby changing the flowfield around the airfoil. 

In computing this flowfield for the propeller case, the 
rotational and induced components of velocity, as well as the 
forward component, must be taken into account as exhibited in 
Figure 2, The reshaping of the airfoil surface by the ice 
accretion and the rough peaks and surfaces commonly present in 
ice accretions, especially in the glaze ice case, serve to 
induce premature transition and often separation of the flow 
around the airfoil, thus spoiling its designed aerodynamic 
characteristics. The associated ice roughness also acts as an 
energy loss device. For the propeller configuration dealt 
with in this thesis, these effects may occur along the entire 
span of the propeller and have the same effect on propeller 
aerodynamics as they are found to have with respect to fixed 
wing aircraft. Ice accretion on a propeller may then 
translate for example into an increase in power required to 
maintain a given flight condition when ice is allowed to form 
on the propeller blades. As a pilot, passenger, or other 
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concerned party it would therefore be desirable before 
undergoing an icing encounter to have some idea of just how 
much power would be required to maintain a flight condition 
with ice accreting on the propeller and other aircraft 
components, so that some decision could be confidently made 
regarding the safety of undergoing or avoiding such an 
encounter. In addition, information provided by this code may 
be used in establishing propeller ice protection system 
specifications, as well as potentially in future icing 
certification efforts. 



CHAPTER 2 


SUMMARY OF PREVIOUS WORK 

Much work has already been done in attempting to better 
.understand ice accretion and ita effect on aircraft 
performance. However, due to the very complicated physics 
which govern the ice accretion process, many questions remain 
unanswered in this area and it is still difficult with present 
technology to accurately predict the degree of performance 
degradation which will result from an arbitrary icing 
encounter. Several recent developments have made the 
prediction of performance degradation, especially in terms of 
drag increase, more reliable and have expanded the range of 
conditions over which relatively accurate predictions can now 
be made. Yet much of the groundwork for these developments 
was laid thirty to forty years ago. This initial work 
involved both experimental and analytical investigations of 
ice accretion, and much of the icing data used in current 
correlation development work was actually collected during 
this early period, from roughly the late 1940's to the late 
1950 ' a . 

A. Early Work 

1. Experimental Effort 

Early propeller performance degradation data was 
gathered using both simulated ice and actual natural icing 
encounters. Corson and Maynard <1> in 1946 used simulated ice 
on propeller blades and measured average efficiency losses on 
the order of three percent, with maximum losses of fifteen 
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percent noted. In 1948, Preaton and Blackman (2) undertook a 
flight test program in which they also noted average decreases 
of roughly ten percent in propeller efficiency due to ice 
formation, as well as the associated increase in drag. This 
study was of special significance due to the fact that an 
attempt was made to examine the effects of ice accretion on 
all major components of the aircraft. Possibly the most 
complete and informative work relative to ice accretion on 
propellers was performed and documented by Neel and Bright (3) 
in 1950. In a flight test program in which efficiency loss 
was measured during actual natural icing encounters, they 
observed losses of roughly ten percent in most cases, with 
maximum losses on the order of twenty percent. The 
performance degradation data of Neel and Bright is still used 
by present researchers to develop or verify analytical 
propeller performance degradation codes. 

Also of note with respect to experimental investigations 
of airfoil icing are a series of test programs undertaken in 
the Icing Research Tunnel at the NASA Lewis Research Center 
and documented by a variety of NACA personnel in the early and 
mid-1950's. Gray and von Glahn <4, 5) looked at the effects of 
ice on the performance of NACA 65A004 and 65^-212 airfoils, 

and leaser increases were noted for both airfoils with rime 
ice accretions. Brun, at. al. (6,7,8) and Brun and Vogt (9) 
produced a series of reports dealing with experimental 
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measurements of performance degradation of NACA 65A004 , 65^- 
'208, and 85^-212 airfoils due to ice formation. These reports 

focused on the impingement characteristics of these airfoils, 
which in turn relate directly to drag increment and decreased 
performance. Included in their investigations were the 
effects of airfoil thickness and angle of attack on droplet 

l • : ■ . : ‘ ■ . ■ . 

impingement values, noting that thickness tends to increase 
the volume of water collected and decrease the rearward limit 
of droplet impingement, and as expected total impingement 
increases with increasing angle of attack. 

Gelder, Smyers, and von Glahn (10) in 1956 investigated 
droplet impingement properties on several airfoil sections of 
various thickness and concluded in part that the total and 
maximum local collection efficiencies were strong functions of 
first of all the modified inertia parameter, a parameter which 
takes into account factors such as freestream velocity, 
droplet density and diameter, airfoil chord, absolute air 
viscosity and air density, and also of thickness ratio and 
angle of attack. 

These early experimental programs not only provided 
valuable insight into the process and effects of ice 
accretion, but also established a data base from which 
analytical ice growth and performance degradation predictions 
could be developed or verified. 

2. Analytical Effort 

In the report by Neel and Bright (3), the authors 
attempted to predict analytically the degree of propeller 
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performance reduction to be expected as a result of ice 
accretion. Using blade element theory as the basis for their 
analysis, and by relating the change in airfoil drag/lift 
ratio to efficiency loss, they predicted efficiency losses 
which agreed at least to the same order of magnitude with the 
experimentally obtained values. 

Bergrun (11) in 1951 offered a method for determining 
droplet impingement characteristics on an airfoil by solving a 
set of simultaneous differential equations which described the 
particle dynamics of a water droplet moving in an air stream. 
This same principle forms the basis for present droplet 
trajectory codes also. Along with Lewis in 1952 (12), Bergrun 
also investigated using probability analysis the atmospheric 
factors responsible for ice formation and identified with a 
quite limited data base the three parameters still felt to be 
of primary importance in the ice accretion process: cloud 

liquid water content, water droplet size, and ambient 
temperature. 

In 1958, Gray undertook an icing study of the NACA 
65A004 airfoil (13) and attempted to develop a drag 
coefficient correlation which would relate ice nature and 
droplet impingement rates to the associated aerodynamic 
penalties for this airfoil. Using experimental icing data, he 
obtained a dimensional correlation which he found to be 
accurate at all but angles of attack greater than four 
degrees, presumably due to the flow separation occurring at 
the higher angles of attack. Gray then published in 1964 a 
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report <14> in which he examined available aerodynamic icing 
data for other airfoils and modified his correlation to make 
it applicable to these other airfoils as well. He 
accomplished this through the introduction of the factor 'r' 
which represents the airfoil leading edge radius of curvature 
in percent of chord. This modified correlation represented 
the only available icing drag coefficient correlation at that 
time, and saw widespread use for several years. Only recently 
has its validity as a general icing correlation been 
questioned, and it will most probably remain in use until a 
suitable replacement ia offered. 

The formulation of Gray's drag coefficient correlation 
represented the last major development with respect to ice 
accretion effects for several years. The advent and use of 
the jet engine and turbofan on the propulsion scene then 
deemphasized the use of the propeller as a propulsive device 
and subsequently drew much attention away from the problem of 
ice accretion, which was of lesser importance for the large 
commercial transports and military aircraft of that time. 
However, with recent increased fuel costs and a renewed 
interest in commuter aircraft and all-weather helicopters, new 
interest in ice protection systems and a need to better 
understand the phenomenon of ice accretion have again come 
about. 
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B. Recent Work 


1. Experimental Effort 

Researchers involved in the renewed study of ice 
accretion recognized the need to expand the quite limited 
experimental data base, especially to include tests of ice 
accretion on the newer airfoils designed for use on current 
general aviation aircraft, propellers, and helicopter rotors. 
Shaw, Sotos, and Solano in a 1982 report (15) discussed the 
results of one such program in which aerodynamic performance 
degradation data was obtained for a NACA 632~A415 airfoil. 

Glaze and rime ice formations were studied in both cruise and 
climb conf igurationa and significant performance degradation 
was noted for all cases. Also evaluated was the effect of aft 
frost growth on airfoil performance, and this phonomenon was 
found to significantly increase section drag coefficient. 

This project was also of importance because it provided new 
data with which to test Gray's drag coefficient correlation. 
The authors found that the correlation was as accurate in most 
cases for the new data as it was for the old data upon which 
the correlation was developed, but it became poorer in 
accuracy for higher liquid water content values. This 
conclusion again pointed out the need for a better, more 
accurate drag coefficient correlation. 

Also in 1982 at the Ohio State University, Bragg and 
Gregorek (16), along with Zaguli (17), used simulated ice 
shapes to investigate various aspects of the ice accretion 
problem. In one test a simulated rime ice shape was applied. 


10 



with and without surface roughness, to a NACA 65A413 airfoil 
and the ensuing performance degradation was measured. Surface 
roughness was determined to be an important factor in 
modelling ice shapes, affecting drag coefficient and C, 

max 

In the second test, they also demonstrated the feasibility of 
using wood shapes to model ice accretions on a NACA 63 fl -A415 

airfoil in order to ascertain the effects of the ice shape on 
the aerodynamic flowfield. 

Flemming and Lednicer in 1983 (18) tested a series of 
scale models of helicopter airfoils to investigate the effects 
of artificial ice accretion on airfoils at high speeds. 
Aerodynamic performance degradation of all the airfoils was 
noted, with the authors citing drag coefficient increases of 
up to roughly 300* in some cases. Also accomplished in this 

l ' i ■ 

teat were better definitions of the boundaries for ice growth 
and ice type in terms of static temperature, Mach number, and 
liquid water content. 

Motivated by a desire to find a relatively simple, 
economical means of collecting performance degradation data 
relative to rotating systems, Korkan, Cross, and Cornell 
(19,20) and Korkan, Cross, and Miller, (21,22) in 1984 
undertook a test program utilizing a model helicopter with a 
53.375-inch diameter, 2.5-inch chord NACA 0012 main rotor with 
a simulated ice shape attached and measured main rotor 
performance degradation. Trends identical to those previously 
seen on larger or full-size airfoils and rotors were noted. 
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with increases of up to 300* in torque coefficient found to be 
required to maintain a given thrust coefficient after the 
addition of the simulated ice to the rotor blades. The 
sensitivity of the tip region of the rotor to the adverse 
effects of ice formation was also seen, with an increase in 
torque coefficient required to maintain a given thrust 
coefficient of as much as 150* measured when the simulated ice 
was extended from the 85* rotor radial location out to the 
100*' location. 

Just prior to this model helicopter test, the issue of 
Reynolds number effects on aerodynamic data was also addressed 
by the same authors (20,22). In a 2-D airfoil test using a 
NACA 0012 airfoil section with a simulated ice formation 
geometrically identical to that used in the model helicopter 
test, the authors collected aerodynamic data for the airfoil 
in both clean and iced configurations over a Reynolds number 
range which included the operating range of the model 
helicopter rotor tip. Reynolds numbers effects appeared 
significant only at the highest Reynolds number tested, 

3.4 x 10 6 . 

2. Analytical Effort 

Much has been learned in recent years as a result of 
these experimental programs, and a major contribution of these 
tests has been in the expansion of the icing data base to 
provide more current data with which to develop and evaluate 
analytical methods for prediction of icing-related performance 
degradation. As in the experimental area, many analytical 
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advances have also been recently made. These advances have 
taken the forma of both better methods of predicting droplet 
impingement characteristics as well as better aerodynamic 
performance degradation correlations. 

Bragg et. al. recently developed a computer program to 
calculate water droplet trajectories and thereby determine 
airfoil impingement efficiencies and theoretical ice shapes 
. . (23). Development of this code represented a major step 
toward the goal of analytical prediction of airfoil 
performance degradation due to ice accretion and the code is 
now often used for this purpose. The code itself will be 
discussed in more detail in a later chapter as it is one of 
the major components of the code development which is the 
subject of this report. Along with this droplet trajectory 
code, Bragg and Gregorek have also formulated a drag 
coefficient correlation for the rime ice condition <16>. In 
this correlation, the change in drag due to ice is given as a 
function of several variables related to flight and 
atmospheric conditions, airfoil geometry, and duration of ice 
accretion. This correlation has seen widespread use since its 
development and has only recently been studied in any detail 
to better define its range of applicability. An unpublished 
investigation performed by this author shows the correlation 
to in several cases significantly overpredict drag increment 
for icing data which appears to be rime-oriented, 
i Cansdale and Gent in a 1983 report (24) detailed the 

development of another two-dimensional airfoil icing code. 
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Unlike the Bragg code, the effects of compressibility, kinetic 
heating, and water runback are taken into account in this 
code, thus making it applicable to both rime and glaze ice 
conditions. Designed to be applied to helicopter 
configurations, the code employs a heat balance analysis to 
calculate the kinetic heating and runback effects. The 
authors have reported good agreement between predicted and 
experimentally obtained ice shapes, temperature distributions, 
and icing threshold conditions. 

Flemming and Lednicer (25) have used experimental icing 
data to formulate new rime and glaze ice drag coefficient 
correlations. In addition, this data was also used to 
formulate separate lift and moment coefficient correlations. 
These correlations, recently published in final form (25), 
have shown some promise in preliminary evaluations and are 
discussed in more detail in a later section. 

Using Bragg's droplet trajectory code, Korkan, Dadone, 
and Shaw (26) have developed a method for predicting 
performance degradation of rotating systems under the 
influence of ice accretions. Limited currently to the rime 
ice condition, the effects of ice formation on helicopter 
rotor or propeller thrust coefficient, torque coefficient, and 
efficiency may be calculated. Originally restricted to the 
helicopter hover mode, extension of the method to include 
forward flight calculations has been performed and documented 
in a recent AIAA paper by Korkan, Dadone, and Shaw (27) . The 
method essentially involves first obtaining non-iced values of 
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angle of attack and Mach number as a function of radial 
location, then determining impingement efficiency and 
accumulation parameter for a given set of icing conditions for 
each radial location, and finally determining the resulting 
drag increment at each radial location. An existing propeller 
or helicopter performance code with the iced drag increment 
input then will produce values of iced performance which may 
be compared with. known or calculated clean values to ascertain 
the. degree of performance degradation which results from the 
given icing encounter. Good agreement between theory and 
experiment has been obtained for both the propeller and 
.helicopter rotor configurations. 

Korkan et . al. have also developed a method of averaging 
, (27) to compute rotor or propeller disk performance in which 
angle of attack and Mach number are first averaged around the 
, disk. From these averaged values, collection efficiency, 
accumulation parameter, and finally drag increment are 
calculated for given icing conditions. This procedure has 
been shown to be virtually as accurate as the previously 
described method but provides a significant reduction in the 
number of necessary calculations and hence computer run time. 

Miller, Korkan, and Shaw in 1983 (28) attempted to 
develop a drag coefficient correlation for the glaze ice 
condition using statistical analysis of icing data as a basis 
of the correlation development. Although a resulting 
correlation was presented, it represented only the state of 
the study at the time and showed a need for much further 
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modification and analysis. Feasibility of using such an 
analysis as the basis for correlation development was 
demonstrated, but the results of such an analysis must be 
interpreted correctly and in tandem with physical observations 
of. the associated phenomenon to make the method truly 
beneficial. 

Bragg (29) has also recently documented a study in which 
he predicted airfoil performance with both simulated rime and 
glaze ice shapes. Various codes were examined for possible 
use in calculating pressure distributions on an airfoil with 
glaze ice attached. The potential flow code of Bristow was 
determined to be suitable for this purpose and provided good 
agreement between predicted and experimental pressure 
distributions. Bragg also investigated the effects of surface 
roughness in the laminar boundary layer on drag increment and 
offered promising results for predicting drag increase by this 
method . 
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CHAPTER 3 


CODES EMPLOYED 

A. Theodoraen Transformation 
The first step in the performance degradation 
determination process is computation of the flowfield about 
the airfoil, section being investigated. The method used in 
this code is a Theodorsen technique of conformal mapping which 
translates the airfoil coordinates to a circle plane and from 
that plane determines the velocity distribution about the 
airfoil, The calculations involved in this method make the 
assumption of potential flow, in which no viscous effects are 
considered. This has been found to be an acceptable 
assumption, allowing the method to yield accurate results, due 
to the fact that viscous effects existing in the boundary 
layer near the leading edge of the airfoil have little effect 
^ on ..droplet trajectory for all but the smallest of droplet 

sizes. The magnitude of the inertia force associated with the 
larger droplets greatly predominates over any viscous 
deflection or interaction phenomena in this region. 

This particular flowfield calculation technique is 
limited to incompressible flow, which is acceptible since 
compressibility effects on droplet trajectory have been found 
to be negligible up to the airfoil section critical Mach 
number for all but the smallest of drop sizes. In recent 
joint studies between NASA and the Royal Aircraft 
Establishment (RAE) <30), back-to-back comparisons were made 
between the impingement efficiency values generated using the 
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incompressible Theodoraen method at NASA and those calculated 
using the compressible flowfield code of Garabedian and Korn 
for supercritical airfoil sections at the RAE. Good agreement 
was observed between the two methods, thereby deemphaaizing 
any compressibility effects and validating the use of the 
Theodorsen transformation in computing the flowfield about 
airfoil sections exposed to icing. 

Theodorsen 'a method makes use of the Karman-Tref ftz 
transformation to relate the airfoil geometry to a near circle 
plane. This particular transformation is used because it has 
been found to yield a better near circle for airfoils with 
finite trailing edge angles by avoiding any cusp atthe 
trailing edge. It is given by the formula 

z = nb <z ♦b) n ♦ <z -b> n 
1 1 

Cz 1 + b) n - <z x -b) n <1> 

where z is a vector which may define some profile such as an 

airfoil section, and z^ is a vector from the origin, 0, to 

some arbitrary point on the circle. Here, b represents a 
vector directed along. the positive x-axis with length OB', 
where B' is a point at which the circle and x-axis intersect. 
The value of n is typically just less than two, and it may be 
observed that for n=2, the Karman-Tref ftz transformation 
simplifies to the more familiar Joukowski profile, given by 



Finally, an iterative procedure calculates the Fourier 
coefficients in the Theodorsen mapping to the exact circle. 
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Using a Newton root-finding scheme, the point in the circle 
plane corresponding to a particular point on the airfoil may 
be determined and from this the flow velocity at that point 
may be obtained given the section angle of attack. 

B. Droplet Trajectory Code 

The droplet trajectory code employed in this analysis 
was developed by Bragg <31> in I960. Using basic theory and 
equations developed several years earlier as its foundation, 
the code extends the earlier work to include the results of 
more recent studies as well as much improved computer 
technology. Essentially, the code is used to compute single 
droplet trajectories and points of impingement on a given 
airfoil, and from these calculations it is able to determine 
doth- total and maximum local collection efficiencies as well 
as a predicted ice shape if desired. 

•U Again, the initial step in the trajectory calculation 
process requires a computation of the potential flowfield 
about the airfoil under investigation. This flowfield 
information, calculated using the Theodorsen transformation 
method described in the previous section, is passed to the 
trajectory calculation portion of the code in the form of the 
aforementioned Fourier coefficients, which yield the flow 
velocity at any given point in the airfoil flowfield for a 
specific angle of attack. With this flowfield information 
available, droplet trajectories may then be computed. 

The trajectory calculation requires solving a 
differential equation which results from the application of 



Newton's Second Law to a water droplet moving in an air 

stream. This equation is given as 

C_ R - - v . 

K T) = D <u-D) + 1 g 

24 Fr 8 g (3) 

where K and Fr are nondimensional parameters known 

respectively as the inertia parameter and Froude number. They 

are given by 

K = 1_&!_u 

18cU (4) 

and 

Fr = U 



From Equation 3 a set of four simultaneous first order 
ordinary differential equations result, which when solved 
yield the droplet trajectory. 

Once the droplet trajectories are known, impingement 
efficiencies may then be calculated. Airfoil local 
impingement efficiency, &, physically represents a ratio of 
airfoil surface droplet mass flux to freestream mass flux, and 
is 



ds <6> 

Here, y fl is the initial particle y-coordinate and S is 

the airfoil surface arc length from the leading edge to the 
point of impact. Figure 3 illustrates graphically the 
calculation of and Figure 4 shows a typical &-S curve. If 
yo is plotted versus S, it may be seen that & is simply the 
slope of the yo -S curve at a given point. Total, or overall. 
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collection efficiency of an airfoil is also an important icing 
parameter and is defined by 



<7) 


where .Yj is /the distance between the initial y-coordinate of 

the upper and lower droplet tangent trajectories and h 
represents the maximum airfoil projected height. Total 
collection efficiency is then a ratio of the amount of water 
‘mass* .collected by the airfoil to the amount in the freeatream 
sector which is swept out by the airfoil. Figure 5 exhibits 
typical variation of collection efficiencies with radial 
,, *.v : /location, .and Figure 6 is an example of the radial variation 
of local angle of attack and Mach number. 

Other capabilities exist within the droplet trajectory 
code , perhaps most significantly the ability to develop a 
theoretical rime ice shape using a time-stepping routine, but 
this function is not necessary in the present study and will 
not be employed or further discussed here. Although it is 
true that the formation of an ice shape on an airfoil will 
indeed alter the flowfield about the airfoil, and hence the P 
curve for that airfoil, this phenomenon is not considered in 
this analysis because the values of & and E used in later 
calculations are taken for the clean airfoil. The importance 
of the collection efficiencies will be manifested in the 
following section dealing with aerodynamic coefficient 
correlations. 
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C. Aerodynamic Coefficient Correlations 
Prior to the initiation of the present study, an effort 
was made to identify all of the dimensionless parameters 
relevant to an icing encounter using Buckingham's Pi Theorem, 
with the goal of identifying and/or justifying the use of 
certain variables in empirical correlations. The Pi theom of 
Buckingham, presented first in 1914, provides a means of 
creating a dimensionless variable set from an Initial group of 
dimensional variables. The theorem states that any function 
of N variables, of the form 


f »Pj »Pj » • ■ .Pfj 5 " 0 

may also be given in terms of (N-K) dimensionless Pi products, 
such that 


f <Pi t ,Pi 2 , . . .Pi NK > = 0 

Here K represents the number of dimensionally independent 
quantities necessary to fully express the dimensions of the N 
variables, and each Pi term is a dimensionless combination of 
an arbitrarily chosen set of K dimensionally independent 
variables. The Pi theorem will be used here to 
nondimensional ize parameters relevant to an airfoil icing 
encounter. 

The procedure is initiated by first listing all of the 
suspected relevant dimensional parameters. In an icing 
encounter, these are taken to be as given in Table 1. 

Table 1. Dimensional parameters considered in the ice 
accretion process. 
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Symbol 

Denotes 

Units 

T 

temperature at which ice 

t 


is formed 


T o 

reference temperature 

t 

: T 

icing time 

T 

w 

liquid water content 

ml” 1 

d 

droplet diameter 

L 

v o 

freestream velocity 

lt” 1 

^ice 

ice density 

ML” 5 

^ a i 

air density 

ML-” 5 


air 


: W absolute air viscosity ML 1 T 1 

v' ; ; c airfoil chord L 

-~i • i .d < t= temperature, T=time, M=mass, L= length) 

Here, ten dimensional variablee (N) and four fundamental 
dimensions <K) are given. Therefore, N-K = 10-4 - 6 Pi 
products will result- It should be noted that although other 
parameters such as total and local collection efficiencies are 
also thought to be important parameters in an icing encounter, 
they are not included in the above variable set because they 
are already dimensionless and would not be affected by the Pi 
theorem which deals only with dimensional quantities. They 
can, however, be added to the final list of nondimensional 
variables which will result from the use of the Pi theorem. 

Certain other physical characteristics of the airfoil 
which may affect ice accretion and resulting aerodynamic 
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performance losses, namely leading edge radius of curvature, 
airfoil thickness, and camber, are not included in the above 
variable set even though they are dimensional quantities. 

This is because the only seemingly meaningful way to 
nondimensional ize these terms is with respect to airfoil 
chord, so that for the purposes of this analysis it will be 
assumed that these parameters are all initially 

nondimensional ized by chord and hence are not Included in the 
variable set. Instead they may be added to the list of 
nondimensional variables at the conclusion of the application 
of the theorem if desired. 

It may be noted that the chord term makes its influence 
felt in the nondimensionalization of these airfoil geometry 
terms, and has also been included in the dimensional variable 
set' to allow for the possible development of dimensionless 
terms which will be consistent with or similar to 
dimensionless quantities already in existence. For example: 

(Accumulation Parameter) 

(Modified Inertia Parameter) 


where 


AC = V 0 UT 
P i ce C 


K_ = f < K , R ) 
0 ' u 


K = P ice d V 0 
18 cU 


R = Vair d 

u 

U 


The next step in the procedure is to formulate a K set 
of four (equal to the number of fundamental dimensions) 
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dimensionally independent variables, from which the Pi 
products .will be developed. To form the set, any four 
variables may be chosen from the dimensional variable list, 
with the restriction that they may not themselves combine to 
form a dimensionless combination. A K set consisting of 
K = <T, T, w, Vo > 

is initially chosen. With this K set, the six resulting Pi 
products are: 


Pi l 

- 

f <T, 

T, 

w» 

Vo , 

To > 

Pi a 

= 

f <T, 

T, 

w. 

Vo , 

d> 

Pi 3 

= 

f (T, 

T, 

w. 

Vo , 

^ice 

Pi , 

= 

f <T, 

T, 

W, 

Vo , 

**air 

P1 B 

= 

f <T, 

T ' 

w * 

Vo , 

u) 

Pi 4 

= 

f (T, 

T, 

W, 

Vo , 

c) 


To create nondimensional terms from these sets, the 
following procedure, as illustrated by an example using Pi t , 

is employed. The Pi product is first written as a product of 
all of its terms, each raised to some power. 

Pi 4 = CT) <T> Q <w) b (Vo) C (To ) d 

It is then rewritten in terms of the fundamental dimensions 
involved: 

Pi 4 = (t) <T) a (ML - 3 ) b (LT - 1 ) C (t) d 
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For this product to be nondimensional ized , all of these ' 
dimensions must go to zero, i.e., the exponents on each ' 
dimension must sum to zero. 

t: 1+d = 0 

T : a-c = 0 

m: b = o 

L: -3b+c = 0 

This set of equations is then solved for a, b, c , and d. 

d=-l b=0 c = 0 a = 0 

The Pi product is then rewritten with these values included. 

Pi t = <T) (To ) 1 = ( T ) . 

To 

which is the first dimensionless quantity. 

Likewise, the other Pi products are found to be 



c 


The resulting set of dimensionless variables as 
determined by the Pi theorem is then 

f < t /To , VoT/d, w/ P , , wVo 8 /U, VoT/c) = 0, 

ice 
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or in the case of predicting drag coefficient change due to 
these variables, 

f (T/To ,- .... . .. ., VoT/c) = AC 

d 

At this point, any other dimensionless variables such as 
collection efficiencies or airfoil geometry terms may be 
included in the equation of desired. The user may then apply 
some correlation development technique to determine the exact 
relationship between these dimensionless variables and the 
associated dependent variable(s) . 

It may be demonstrated that the resultant set of 
dimensionless variables is dependent on both the set of 
dimensional variables which is used and on the K set which ia 
initially chosen. Since there are a finite number of 
dimensional variables and a finite number of K sets which can 
"be chosen, there will be a finite number of dimensionless 
variables which will result for a given problem. However, not 
all of these variables will necessarily be produced by one K 
set. For the icing encounter problem discussed here, this can 
be illustrated by choosing a second K set: 

K * < T, w, c, Vo ) 

This particular K set is especially noteworthy for the icing 
encounter problem since it contains parameters which are 
thought to be of primary importance in the ice accretion 
process: temperature, liquid water content, freestreah 

velocity, and airfoil chord length. The Pi products for this 
set are: 
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Pi i 

S 

f <T, 

w. 

c. 

Vo , 

To > 

P1 a 

= 

f <T, 

w. 

c. 

Vo , 

T) 

p 1 J 

= 

f <T, 

w. 

C, 

Vo , 

d> 

P1 H 

= 

f <T, 

w. 

C, 

Vo , 

^ice 

Pi « 

= 

f ct; 

w. 

C, 

Vo , 

P air 

Pi * 

s 

f C T , 

w. 

C, 

Vo , 

U> ' 


Following the procedure previously outlined for 
nondimensionl izing these products, the resultant set of 
dimensionless variables is: 

f<T/To, Vo T/c, c/d, w/ P , w/ P , wcVo/W) = 0, 

ice air 

which differs somewhat from the set produced by the other K 
set used. So for any given problem it is necessary to 
consider all of the possible K sets to obtain a complete list 
of dimensionless variables which may be produced. Knowledge 
of the physics of the problem under consideration should then 
provide a starting point in the analysis of these variables as 
they relate to the given dependent variableCs). 

In the icing encounter problem, it was desired to use 
the Pi theorem to provide some justification for the use of 
the terms AC and Ko , previously defined, in the development of 
a glaze ice drag coefficient correlation. Although the Pi 
theorem does not produce these terms directly, it does provide 
the following dimensionless terms: 
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< 1 ) Vo T/c 

< 2 ) w/P. 

ice 

(3) PiceV Q d 

M. 


< 4 ) d/c 

'(,5.) PairV 0 d 
W 


' It may be noted that AC is simply the product of terms 
<1) and (2). Ko is a function of K, the product of terms (3) 
and (4), and R , which is identical to term (5), e.g., 

V 

AC = OT w 


< 1 > 


< 2 ) 


K = 


p , v -d 

ice 0 
P. 

(3) 


d 

c 


P, V o d 

ice 0 

W-c 


(4) 


R = PairV d 
u Q 

• ' 1 U 

<5) 

Should the reader wish to verify that these terms are indeed 
obtainable from the original set of dimensional variables, 
they can be found as follows: 

(1) - from K set (T, T, w, Vo) 

(2) - from K set <T, T, w, Vo ) 

(3) - from K set (T, d, Vo , P , ) 

ice 

<4) - from K set (T, w, c, Vo > 

(5) - from K set <T, d, Vo , P . > 

air 
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It can be aeen that the Pi theorem does indeed provide 
some justification for the use of the variables AC and Ko in 
developing a correlation, and any correlation which does 
result will be in terms of dimensionless variables obtainable 
by Buckingham's Pi theorem. 

At present, the only means available for analytically 
determining performance degradation of an airfoil in terms of 
its aerodynamic coefficients is through the use of empirical 
correlations. As previously indicated, several individuals 
have attempted development of one or more of these 
correlations and generally have met with only limited success 
in terms of accuracy and range of application. Incorporated 
into the present study are the two correlations currently most 
widely used to predict drag increment, as well as a aeries of 
new correlations. These empirical correlations represent only 
an interim solution to the problem of assessing aerodynamic 
performance losses due to ice formation. A potentially more 
precise solution technique currently under development 
involves a Navier-Stokes solution for the flowfield around an 
iced airfoil <320 . 

1 . Gray Correlation 

The first correlation, formulated by Gray (14) roughly 
twenty years ago is of the form 

AC= [8.7 x 10~ 8 Tu W6 (32-To)°’ J ] (1+6 { 1 + 2 . 52r° ' 1 sin* 12oO 

D <±— max 
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) 3 


sin 2 C543 w ( E )i /J - 81 +65.3 ( 1 - 1 

32-To 1.35* 1 1 + 35°* 

-0.17 sin 4 llctf) ) (8) 

r 

where here AC^ represents the actual numerical value of drag 

coefficient increment; i.e., AC. = C,(iced>- C, (clean) . As 

d d d 

can be seen, this correlation involves combinations of various 
dimensional parameters, and the dimensionless total and local 
■collection efficiencies calculated in the droplet trajectory 
code previously discussed also appear. This correlation also 
has the capability of computing drag coefficient increment at 
angles of attack other than that at which the ice was formed, 
through the use of the c* and en^ terms. Here c/. represents the 

angle of attack at which the experimental data was taken, or 
the angle at which the effects of the ice accretion are to be 
evaluated, and <x^ represents the angle of attack at which the 

ice was formed. 

The dimensional variables have not in all cases been 
grouped so as to form dimensionless quantities, and this has 
evoked some criticism by various persons since the 
correlation's development. This criticism is based on the 
idea that any correlation which involves dimensional 
quantities will be to some extent a function of the data set 
from which it was formulated. Thus, if the correlation was 
applied to some other set of data in which the range of one or 
more of the relevent parameters differed from that of the 
original data set, the variable affected may have an undue 
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effect on the prediction which was not accounted for in the 
original, limited data set. Indeed, as is seen in Figure 7, 
Gray's correlation has been found to predict quite poorly for 
much of the recently generated NASA glaze ice data, which is 
typefied by significantly higher liquid water content values 
than are found in the data set used by Gray in developing his 
correlation. Nevertheleaa, since no other glaze ice drag 
coefficient existed until quite: recently, Gray's correlation 
still sees widespread use and has been included in this 
program as an option available to the user. 

Because Gray's correlation accounts for only changes in 
drag coefficient, modification of has been added to take 

into account the decambering effect of ice accretion, 
resulting in a shift in the -<* curve and a decrease in 

max 

due to flow separation on the airfoil upper surface. Trial 
values tested by the author have indicated that reducing to 

95* of its clean configuration value yields an acceptable 
approximation for iced . This factor is incorporated into 

the present study for both Gray's and Bragg's correlations. 

2. Bragg Correlation 

A correlation has been developed by Bragg (16) for the 

rime ice condition. It has the form 

AC . = 0.01 <15.8 ln<k/c> 28000 A E * I ) <9> 

d c 

where Ac is accumulation parameter, given by 

Ac = uwT 

P, c <10> 

ice 
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and I represents a drag constant which variea according to 
airfoil type as shown in Table 2. 

Table 2. Constants in Bragg Drag Coefficient 
Correlation 


Airfoil Type 
NACA 4 and 5 digit 
NACA 63 series 
NACA 64 series 
NACA 65 aeries 
NACA 66 series 
In this correlation, 


Drag Constant. I 
184 
218 
232 
252 
290 

AC ' represents a fractional change 
d 


in, rather than an actual numerical value of, drag 
coefficient. Here AC d may be used to find the iced airfoil 

drag coefficient using 

C, = (ltACJC^ (11) 

d. , d d , 

iced clean 

Previous studies of this correlation applied to NACA and 
NASA 2-D icing data have indicated a strong tendency for the 
correlation to overpredict drag increment. Because the 
predictions are commonly off by a relatively constant 
percentage with respect to the Neel and Bright propeller 
performance data, a factor was introduced which reduces the 
degree of drag coefficient increase calculated by the Bragg 
correlation. Replacement of the initial constant in the 
equation, 0.01, by 0.0008 produced much more acceptable 
predictions and was therefore incorporated into the present 
analysis. 
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/ 


Bragg has also recently proposed a new, modified form of 
his correlation (33) which has also been Included in the 
present study. This new form of the correlation given below, 
appears to work well for the cases examined to date, but it 
has not yet been applied to the 2-D NACA and NASA icing data 
with which the other correlations were tested. 

AC . = 0.01 ( 15 . 81n ( k/c ) ♦ 1171 AcE ♦ I) (12) 

d 

3. Flemming Correlations 

The third correlation, or more precisely set of 
correlations, are those developed by Flemming (25). Whereas 
Gray and Bragg dealt strictly with drag coefficient 
correlations, Flemming has also formulated correlations for 
lift and moment coefficients for both glaze and rime 
encounters. These correlations are of the forms: 

For a glaze ice encounter, 

AC, = KD. [.00686 Ko < t/c) 1 * S (<*+6) - . 0313 ( r/c) 2 -»-KD ( . 006 ) M 2 ' H 3 

a * 

( w Tc.( c )°* 2 /( c 3 1 ' 2 3 C1-8AC ( V helo)] (13) 

^ 

.1524 .1524 278 

For a rime ice encounter. 


AC . = [ . 1581n(k/c) -*-175( V )WT E ♦ 1.703 ( <X+6 ) 

. d c 

P. c 10 

ice 


1 1 -8AC ' < help ) 3 C, 

278 clean 


and for both glaze and rime encounters, 

AC. = (-.01335Ko (t/c) Co< + 2 + KL 4 ( .00555) (<x-6) 8 3KL3 CWT 
% 1 c 

( c ~ ) 0 * 2 / ( c ' ) 1 * 2 3 (15) 

.1524 .1524 
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and . 


ACm = [< .00179-. 0045M) < . 00544 > Ko */ < t/c> 8 * 7 

.00383M< 1-63. 29^0 3 (16) 

[WT < c )° * 8 /< c 3 1 * 8 3 
c 

.1524 .1524 

where Ko is the modified -Inertia parameter, given by 

Ko = K < 1 ■*■ . 0967 R - A597 > (17) 

u 

and KL, KL 4 , and KD, and KD 4 are functions of temperature and 

angle of attack and are defined within the subroutine. 

These correlations were developed using data obtained in 
a series of tests conducted in the Canadian National Research 
Council's High Speed Icing Wind Tunnel in late 1982- The ten 
airfoil scale models tested, with typical chorda of roughly 
six inches, were designed primarily for use in helicopter 
applications and cover a wide range of helicopter airfoil 
shapes. For correlation development purposes this range 
should enable the resulting correlations to adequately account 
for airfoil geometry as a factor in aerodynamic coefficient 
increments due to icing. 

The correlations are, however, dimensional in nature and 
again were developed using model data for high speed 
applications. The full effects of these two items on the 
range of applicability of the correlations remains to be seen. 
The correlations have been found in a recent study to 
underpredict AC^ in many cases however. In spite of this. 
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Flemming's correlations still hold promise and do in fact work 
well in many cases. 

4. Miller Correlation 

A fourth general drag coefficient correlation 
development has been attempted by Miller, et. al. Described 
in detail in Reference 28, this effort used as the basis for 
the correlation formulation the method of statistical analysis 
as applied to a set of experimentally obtained icing data. 
Although a correlation was presented in Reference 28, it 
represented only a status report and not a final result. 
Further study and modification of this correlation is 
necessar,y before it may be deemed useable, and no 
modifications made to date have proven satisfactory. The 
method of statistical analysis is a powerful tool in 
correlation in development and should yield greater success at 
some future date, but because the present form of the 
correlation is not acceptable in terms of accuracy of 
prediction it will not be incorporated into the code. 

So there are presently available three correlations or 
sets of correlations which are integrated into the code. They 
provided values of drag, lift, and/or moment coefficient 
increment for both glaze and rime ice encounters. These 
increments may then be passed to a propeller performance code, 
described in the following section, to evaluate the ensuing 
propeller rotor performance degradation. 
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D. Propeller Performance Code 
The .. propel ler performance code which is integrated into 
the performance degradation program is based upon a linearized 
inflow propeller strip analysis, developed by Cooper (34) in 
1957. Strip analysis involves the calculation of aerodynamic 
forces at selected spanwise blade locations for a. given 
operating condition and propeller geometry. From these 
differential forces, given in terms of thrust and torque 
coefficients by 


riC c 0 

T = k - M (C, cos+-C sin+) 
D x A d 

dx 


(18) 


where k= 900a b 

n“ d 


(19) 


and 


dC 


dC_ 


_J0 = x 

dx 2 dx 


C, cos^+C, sin* 

u > 

C. cos4— C ,sin4 > 

A a 


( 20 ) 


thrust and torque coefficients may be obtained by integrating 
these terms along the blade. Propeller efficiency and power 
absorbed may then be calculated by 


t) = J C T = ( V ) C T 


( 21 ) 


nD 21TC. 


Q 


and 


BHP = C o P * V 


( 22 ) 


550 


Cooper's propeller performance analysis is unique due to 
the fact that it assumes, or actually approximates, a linear 
induced flow distribution on the propeller blade. Typically 
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an iterative procesa ia required to obtain the induced angle 
of attack and' lift coefficient at a apecific radial location 
aince they are interrelated. However, Goldstein (35> haa 
solved for the radial distribution of circulation for a 
lightly loaded propeller having a finite number of bladea and 
has established the relationship between and . Cooper 


has then approximated the 


8 C. 

P * 


curves by straight lines 


and haa plotted the slopes of theae lines versus advance ratio 
for various spanwiae locations. From theae plots, and using 
available aerodynamic data, it ia then possible to determine 
values for an induced angle of attack. When experimental and 
analytical results obtained using this method were compared, 
agreement to within three percent was seen for over ninety 
percent of the cases investigated. 
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CHAPTER 4 


PERFORMANCE DEGRADATION CODE ASSEMBLY 
■ In order to integrate each of the codea involved into a 

single performance degradation code, heretofore known aa 
ICEPERF, several phaaea of modif icationa were required. These 
modif ications were relatively minor in that they did not 
change the baaic functions, operation, or logic of any of the 
codea but rather aerved to make the tranaition from one 
aection to the next more aimple and automatic. Thia 
automation then eliminated the neceaaity of any uaer 
intervention in the run atream. 

The uaer-f riendly nature of the code waa a primary 
objective of thia effort, and although aeveral filea are 
necessary to run the code very few changes are required from 
one run to another, for the aame airfoil. The varioua 
components of the code were selected on the basis of 
availability and accuracy of results. Each component waa run 
and tested separately and extensively to demonstrate its 
capability to consistently function properly, and all 
displayed good performance in these runs. Also of importance 
in the component selection and modification was a desire to 
maintain a high degree of modularity in the code, such that a 
different flowfield computation, new correlation, different 
propeller analysis, etc. could be inserted with a minimum of 
changes to the code. This also permits the code, now applied 
to propellers, to be relatively easily extended to a 
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helicopter configuration. A flowchart denoting the current 
major elements of the program is provided in Appendix 1. 

The first phase in the code integration process involved 
combining the flowfield and droplet trajectory codes 
previously described. Because the droplet trajectory code was 
designed to operate using much of the output from the 
Theodorsen transformation code, this step was relatively 
straightforward. The modified code has been set up to handle 
a maximum of four input propeller radial locations, and 
current runs of the code have been made with four input 
stations, at the 30, 50, 70, and 90 percent radial locations. 
An investigation of the benefits of using more input stations, 
offset by the additional cost due to the extra run time 
required, has not yet been performed. 

Next, a set of subroutines was created to predict drag 
increment due to ice accretion, given the appropriate 
impingement efficiency values output from the trajectory code. 
The correlations of Gray, Bragg, and Flemming are included 

in the subroutines, with Flemming's correlation also being 
capable of determining lift and moment coefficient increments. 
The user may select any of these correlations in each run of 
the program, and the choice of correlation should be dependent 
on the type of ice formation expected for a given input 
condition. 

The correlations of Gray and Bragg use values of total 
and maximum local collection efficiencies in computing drag 
coefficient increment. One value for each of these parameters 
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is generated for each radial location input, ao that a maximum 
of four sets of collection efficiencies are available as 
output from the trajectory calculation section. However, it 
is best when using the propeller performance code to include 
considerably more input radial locations in the calculations 
(eleven stations are used in the analysis of a propeller 
utilizing double-cambered Clark-Y airfoil sections which 
follows in Chapter 5). In order to obtain values for drag 
coefficient increment at each of these locations, the 
collection efficiencies at these stations must be known. 

Rather than running the Theodorsen transformation and droplet 
trajectory codes eleven times to obtain collection efficiency 
values at each of these stations however, a curve fit routine 
is employed which takes the values computed at each of the 
four input radial locations and fits a curve through them. 
Flemming's correlations compute values of collection 
efficiencies empirically, so that it is not necessary to run 
the Theodorsen and droplet trajectory sections of the code if 
Flemming's correlations are to be employed. 

In this code, a cubic spline fit routine is used to fit 
separate curves for total maximum local collection 
efficiencies between the 30* and 90* radial locations, or 
between the two most extreme sections input by the user. This 
method has been found to yield good results due to the fact 
that the collection efficiencies are both smooth, continuous 
functions or propeller radial location. Due to the nature of 
the spline fit routine, it is applicable only between the 
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innermost and outermost radial locations input. Outboard of 

the outermost input station, and inboard of the innermost 

station, some other technique is required to compute values 

for the collection efficiencies as a function of radial 

location. In this analysis a linear fit, computing the slopes 

at the endpoints of the spline fit and extending straight 

lines inboard and outboard of these endpoints, is used to 

approximate the curves where the spline fit is not. applicable. 

With this method, values for both total and maximum local 

collection efficiencies are available for any radial location. 

Cooper's propeller performance code has been utilized to 

calculate values of propeller thrust, torque, and efficiency 

for a given condition. The previously mentioned drag 

increment subroutines have been incorporated into the 

propeller performance code to modify the drag coefficient and 

values used by the code in computing iced propeller thrust and 

torque. In Equations 18 and 20, the C, and C . terms are 

Ad 

modified by the factors (1+AC^) and (l+AC d > respectively, such 
that 

< clean ) * <1 + £0^) <23) 

and 

C =C , (clean) *<1 +AC_.) <24) 

d d d 

so that iced propeller thrust and torque coefficients as well 
as propeller efficiency may then be computed. The code first 
calculates clean, or non-iced, values for these quantities, 
and then calculates using the user-selected aerodynamic 
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coefficient correlation the iced values for these same 
parameters. Various input parameters have been altered or 
omitted from the original code to make it more conducive to 
the purposes of the present program, but the esssence and 
mechanics of the code remain unchanged. 

The entire code then operates by first computing clean 
propeller performance values for one advance ratio. Next, the 
Theodoraen transformation is performed, or previously 
calculated results are read for the first input radial 
location. The droplet trajectory section then computes 
impingement efficiency values for the radial location using 
the flowfield from the Theodoraen transformation and a section 
angle of attack computed in the propeller performance section. 
Impingement efficiencies are then calculated for the remaining 
input radial stations following the same procedure. The code 
spline fits the impingement efficiency curves and computes 
using these values and the selected aerodynamic coefficient 
increment correlation the iced propeller performance data. 

This entire procedure is then repeated for each additional 
advance ratio input. 
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CHAPTER 5 


RESULTS 

In the Neel and Bright teat, a C-46 twin-engine aircraft 
was used to obtain propeller efficiency losses as a result of 
ice formation. Ice was allowed to form on the propeller of 
the right engine while the left engine propeller waa kept 
clean. Following this procedure, data was recorded for twelve 
natural icing encounters and typical efficiency losses of four 
to ten percent were noted. The propellers on the C-46 
utilized double-cambered Clark-Y airfoil sections, for which 
airfoil coordinates used in ICEPERF were obtained from Borst 
(36) . 

The Neel and Bright experimental data provided a data 
base with which to evaluate the ICEPERF iced propeller 
performance code. For ten of the twelve encounters documented 
in the test program report (Table 3), calculations were 
performed for the propeller in both clean and iced 
configurations. Results and comparison of analytical and 
experimental performance values are here presented for each of 
the three correlations employed, specifically those of Bragg, 
Gray, and Flemming. It should be noted that in all 
performance calculations presented in this report, 
consideration was given to the radial extent of icing such 
that no ice was allowed to accrete analytically outboard of 
the experimentally observed radial icing extent. These 
results therefore differ from previous calculations shown by 
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this author <37) in which.no restriction was placed 
analytically on the radial extent of icing. 

The results presented in this paper may also be compared 
directly with calculations performed by Korkan, et . al (27). 
Their work involved generation of propeller performance 
predictions for the Neel and Bright data as well, using the 
Theordorsen transformation flowfield code, Bragg droplet 
trajectory code, and Cooper propeller performance codes 
separately. In this initial study by Korkan, et. al . , 
pressure altitude was input and the standard atmospheric 
pressure, temperature, and density were then used in the 
calculations rather than the experimentally obtained values 
for these parameters. The use of actual flight conditions, as 
employed in all calculations presented in this report, was 
found to alter the computed thrust and power coefficients by 
roughly one percent. Potential differences in blade angle 
specification from study to study may also have a significant 
impact on the predictions. Finally, the initial constant in 

- • t • . . 

the original Bragg correlation was decreased in both studies 
in an attempt to produce reasonable performance values. This 
then introduces additional variation between the corresponding 
iced calculations. 

Clean performance predictions were found to agree well 
with the experimental data for all of the encounters 
investigated. (Tables 4-6 summarize the analytical and 
experimental data for Encounters 2, 7, and 12 respectively.) 

In all cases, calculated thrust and power coefficients were 
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found to agree with the corresponding experimental values 
within an average of 12.5 percent. and for a given 

encounter and advance ratio tended to be underpredicted by 
roughly the same percentage, so that resulting propeller 
efficiencies compared with experimental measurements to within 
six percent. Since propeller efficiency is essentially a 
ratio of thrust coefficient to power coefficient however, 
errors of similar magnitude in C^. and should then cancel so 

that the resulting efficiency values would be much less 
sensitive to error in performance calculations and would not 
be as good an indicator of analytical prediction accuracy. 
Therefore statements herein relate primarily to thrust and 
power coefficient calculations rather than propeller 
efficiency. 

The results shown in this report all involved use of the 
actual blade angles as measured and reported in the 
experimental test program. It was found however that by 
increasing the reported blade angle setting by 0.7 degrees, 
much better agreement between theory and experiment could be 
obtained for all encounters. Figures 8-10 are illustrations 
of the effect of blade angle variation on clean thrust and 
power coefficients for Encounter 2. This anomaly could be an 
indication of a deficiency in the performance calculations, 
but because a constant increase in blade setting produces much 
better correlation with experiment for several encounters over 
a range of flight conditions, it more likely signifies the 
existence of a bias error in the experimental measurement of 
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blade angle. These calculations serve also to indicate the 
sensitivity of propeller performance predictions to propeller 
geometry. 

The code was found to often overpredict the clean 
propeller efficiency at the higher advance ratios.. Typically, 
performance predictions agreed with experiment to within a few 
percent at the lowest advance ratios and gradually worsened to 
the point where they were often off by roughly twenty percent 
at the highest advance ratios investigated. It is postulated 
that this may be due to an inability of the code to adequately 
account for compressibility effects, which become more 
significant as advance ratio increases (via increased forward 
velocity for a fixed rpm). This effect should be studied 
further in subsequent efforts. 

When the restriction on radial icing extent was imposed 
in the present study's calculations, better agreement with the 
experimental data was obtained than was previously presented. 
Currently there is no mechanism in the code for determining 
the actual radial icing extent. Rather, the user must specify 
the desired or assumed radial location inboard of which iced 
calculations are to be performed. In this study, approximate 
radial icing extent for each encounter was reported by Neel 
and Bright and these values were then input for the iced 
performance calculations. The following discussion details 
results obtained by each correlation for the several 
encounters investigated. 
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A . BRAGG 

Two forms of Bragg's basic rime ice correlation are 
available in ICEPERF, as discussed in a previous section. The 
first of these is a modified form of Bragg's original 
equation, in which the leading constant was changed from 0.01 
to 0.0008 in an effort to produce more realistic drag 
coefficient increments. This change was made after initial 
runs of the code indicated that the original equation was 
producing quite unreasonably large performance changes due to 
ice formation, as a result of inordinately large predicted 
AC^ values. The 0.0008 constant was arrived at by best- 
fitting the correlation's performance predictions to Encounter 
2 data and was then found to perform acceptably for the other 
encounters as well. Future use of the correlation may suggest 
additional modifications of this constant. 

f 

This modified form of Bragg's original correlation was 
then applied to all of the available Neel and Bright data 
regardless of assumed rime or glaze ice type, as was Bragg's 
new form of the correlation. Reasonably good agreement 
between experiment and theory was obtained using both 
correlations for all encounters. In all cases the original 
but modified Bragg correlation produced more accurate results 
than did the newer correlation, but the differences were 
generally small. The fact that the newer, unmodified 
correlation gave results almost as accurate as those of the 
older form which was specifically modified to better fit the 
propeller data indicates that the changes made by Bragg in 
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forming the new correlation improved the quality of the 
correlation significantly. Both correlations handled the less 
severe icing encounters (icing extents below 60X) well (i.e.. 
Encounter 7, Figures 14-16) , much better than did the Gray or 
Flemming correlations. In the most severe case however 
(Encounter 12, Figures 17-19) both correlations seriously 
overpredict the degree of degradation, especially in terms of 
C , and produce C -J curves which inflect upward at the lower 

p p 

advance ratios rather than downward as would be expected. 

It should also be noted that Bragg's equations are for 
AC^ only. Lift coefficient will also generally decrease as a 

result of ice formation at the leading edge as the flow there 
is spoiled to some extent, but this is not accounted for in 
the Bragg equations. A factor was therefore introduced to 
provide a simple representation of this phenomenon by 
assigning the iced section lift coefficient a value equal to 
95X of its clean value. Although this seems to be a 
reasonable typical value for AC^ , it would be desirable to 

incorporate a more discriminating form of determining AC^ in 
future calculations. 

B. GRAY 

As was the case with the Bragg correlations, the Gray 
correlation also provides its most accurate performance 
predictions when radial icing extent is small. It is 
difficult to make any inference regarding the quality of the 
correlation under such conditions due to the fact that the 
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propeller Is loaded most heavily at the tip and its 
performance will be moat sensitive to geometry changes (i.e., 
ice accretions) in the tip region. Farther inboard, even 
quite severe ice growths will have much less effect on the 
overall performance and consequently any correlation, 
regardless of how great its predicted A values may be, will 

appear to work well in low radial icing extent cases. Thus 
the true measure of the quality of these correlations comes in 
the more severe icing encounters such as Encounters 2 and 12. 
In these two encounters, as can be seen in Figures 11-13 and 
17-19, the Gray correlation fails to accurately predict the 
iced performance, yielding values which overpredict the degree 
of degradation. Indeed, upon examining the actual AC rf values 

predicted by Gray for Encounter 12 it was found that this 
correlation was predicting drag coefficient increases of 
typically 300* at lower advance ratios to as much as 6000* at 
the outer stations at high advance ratios. The Gray 
correlation, like that of Bragg, also predicts only drag 
coefficient increment, so the same five percent reduction in 
for the iced configuration was again employed when Gray's 

correlation was used. The Gray correlation itself was 
unmodified for use in this study, and it is likely that 
introducing some constant into Gray's equation to affect a 
constant percentage decrease in its calculated AC^ values, as 

was done with the original Bragg equation, would greatly 
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improve its ability to calculate iced drag coefficient 
increment and hence iced propeller performance. 

C. FLEMMING 

The Flemming correlations were applied twice to each 
encounter, once with icing restricted by the user-input value 
and then allowing the Flemming correlations to predict the 
radial icing extent. The second runs were performed to check 
the no-icing capability of the correlations and to ascertain 
whether this capability would properly reflect existence of 
ice under the influence of 3-D, rotating forces (kinetic 
heating and shedding). Instead of properly predicting radial 
icing extent however, this feature dictated no icing in the 
first 30* of propeller radius for any encounter, and then 
overpredicted the radial icing extent outboard of the 30x 
station for eight of the ten encounters (Table 7) . 

In the case of the user-input icing extent, the no-icing 
problem manifests itself in an underprediction of the 
propeller performance degradation for all encounters. As the 
figures illustrate, many of the Flemming-predicted iced thrust 
and power coefficient values are identical or virtually 
identical to the clean values, indicating analytically little 
if any degradation. While these results may appear reasonable 
at first glance in cases where the actual degradation is 
minimal, they actually serve only to mask the underlying 
problem and reason for the negligible degradation prediction. 
It should also be noted that in Figures 11 through 19, in 
cases where the clean theoretical curve is not easily visible. 
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it has been hidden by the Flemming correlation curves which 


are often virtually coincident with the clean data curve. 

When no restriction is imposed upon radial icing extent, 
Flemming's correlations still tend to underpredict the degree 
of iced performance degradation in several cases. This occurs 
even though Flemming's correlations predict more extensive 
icing radially than was seen experimentally, which would 
indicate that the percentage decreases in thrust and power 
coefficients predicted by Flemming are actually smaller than 
the experimental values. For the encounters in which the non- 
restricted Flemming correlations do predict C^, and well. 


the predicted radial icing extent is always greater than the 
measured extent, sometimes by as much as 60%. 

It would therefore appear that the Flemming correlations 
have two significant deficiencies. The first of these is the 
no-icing calculation which has failed in this study to 
properly indicate either ice existence or icing extent. In 
fairness to the correlations it should be pointed out though 
that various forces not stressed in the 2-D correlation 
development (i.e., centrifugal forces leading to ice shedding) 
may have a significant effect on the radial icing extent, and 
in only one out of ten encounters did the Flemming 
correlations actually underpredict the extent of ice. The 
second apparent shortcoming of the correlations is the 
underprediction of section drag and possibly lift coefficient 
increments. This same trend toward underpredict ion of section 
has also been seen in application of the Flemming 
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correlations to available NACA and NASA 2-D icing data in a 
separate study (38) . 
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SUMMARY AND CONCLUSIONS 


Application of the ICEPERF code to the Neel and Bright 
experimental data baae has provided much insight into the 
code's performance prediction capability for both clean and 
iced configurations. The code has been shown to run properly 
and to produce realistic thrust and power coefficient values, 
especially when radial icing extent is known and input. As a 
result of this study several areas of future development and 
potential improvement, discussed in the following paragraphs, 
have been identified. 

Currently restrictions and simplifications exist within 
the code which, if removed or modified, should yield better 
performance evaluations. Development of a correlation to 

be used in conjunction with the Bragg or Gray correlations 
would provide a more precise and complete representation of 
the effects of ice formation on section aerodynamics. 

Secondly, incorporation of some computational mechanism of 
predicting the transition from rime to glaze ice, or vice 
versa, on the blades would enable the code to use the proper 
correlation at each radial computational station, rather than 
assuming all rime or all glaze ice in any given encounter. 
Because of the generally moderate radial icing extent and only 
slight degree of ensuing degradation, this change will 
probably have only a very minor effect on the predicted 
propeller performance. This is in contrast to the helicopter 
rotor icing situation in which radial icing extent and the 
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related performance degradation are typically more severe for 
a given icing condition. 

Another potential area of improvement lies in the 
accuracy of the experimental data used for comparison with 
analyaia. One auch area with reapect to the Neel and Bright 
data is the actual radial icing extent for each encounter. Aa 
apecified in the report, the true extent ia somewhat unclear 
and may possibly vary by aa much aa 20 %. Values of radial 
icing extent presented in this report represent a beat 
estimate of the extent aa interpreted from the information in 
Neel and Bright's report. 

The present study has also provided insight into the 
type and accuracy of experimental data which would be 
desirable in future test programs. Additional propeller icing 
data of any type would be useful in further validating and 
analyzing the code at the present time. Future teats though 
should include, in addition to flight and atmospheric 
parameters and icing time, specification of icing extent as 
accurately as possible. Also desirable would be some 
photographic or other indication of the ice type so that the 
proper correlation could be applied at each radial location. 
Finally, sensitivity of propeller performance to blade angle 
setting also has indicated the need for quite accurate blade 
angle measurements. Certainly it is desirable to have a high 
degree of accuracy in the thrust and power measurements as 
well . 
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With respect to the empirical lift and drag correlations 
used in this code, results obtained in this study concurred 
with previous work done involving these correlations in 
indicating deficiencies which exist in each correlation. 
Overall, the Bragg correlation < s) produced the most reasonable 
iced propeller performance predictions but overprediction of 
drag coefficient increment by Bragg' a correlation (a) remains a 
problem, as illustrated in the more. severe icing encounters. 
Overprediction also remains a problem with Gray's correlation, 
whereas Flemming's correlations typically underpredict the 
aerodynamic coefficient increments. It is evident then that 
much more 2-D validation of all these correlations is 
necessary before iced propeller performance can be predicted 
consistently using this method and before results of this 
analysis can be fully evaluated. It should be realized that 
any ultimate analytical treatment of ice accretion on rotating 
systems will potentially involve much less empiricism than is 
currently required and that the use of these correlations 
represents an interim approach to the problem, which will 
continue to be used until more sophisticated analytical 
procedures such as Navier-Stokes or interactive boundary layer 
analyses of iced airfoil flowfields become available. 
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FIGURE 1. The two basic ice types 
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FIGURE 2. Velocity components of propeller section 
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FIGURE 4. Variation of local impingement efficiency with 
airfoil surface arc length; r/R = 0.9, J = 1.18 
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FIGURE 5. Radial variation of accumulation parameter, maximum 
local collection efficiency, and total collection 
efficiency for Encounter 2, J = 1.18 
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FIGURE 6. Radial variation of local angle of attack and Mach 
no. for Encounter 2, J = 1.18 
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FIGURE 7. Comparison of experimental results with Gray's 
predicted AC, values 
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FIGURE 10. Effect of propeller blade setting on predicted 
propeller efficiency. Encounter 2 
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FIGURE 11. Variation of thrust coefficient with advance ratio, 
clean and iced (iced to r/R = 0.7), Encounter 2 
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FIGURE 12. Variation of power coefficient with advance ratio, 
clean and iced (iced to r/R = 0.7), Encounter 2 
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FIGURE 13. Variation of propeller efficiency with advance ratio, 
clean and iced (iced to r/R = 0.7), Encounter 2 
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FIGURE 15. Variation of power coefficient with advance ratio, 
clean and iced (iced to r/R = 0.4), Encounter 7 
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FIGURE 17. Variation of thrust coefficient with advance ratio, 
clean and iced (iced to r/R = 0.9), Encounter 12 
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FIGURE 18. Variation of power coefficient with advance ratio, 
clean and iced (iced to r/R = 0.9), Encounter 12 
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FIGURE 19. Variation of propeller efficiency with advance ratio, 
clean and iced (iced to r/R = 0.9), Encounter 12 
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Propeller Blade Pressure Avg , Liquid Avg. Droplet 

Encounter Propeller RPM Angle at 0.75R,deg. Altitude ft. Icing Time, min. Air temp.,«F Water Content,?/*" 1 Diameter , microns 
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Table 4. Conparlaon of predicted and experimental propeller performance for 
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Table 7 


. Experimental and predicted (Flemming correlations) radial 
extent for each encounter. 
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APPENDIX 2 


USER'S GUIDE 

ICEPERF has been designed such that user intervention in 
the run stream has been eliminated and the amount of user 
input is minimized. Still, with the wide scope of work and 
calculations which the code must perform, a substantial degree 
of input is required of the user. This information is 
provided through several files which are read at various 
points in the program and which are separated such that each 
contains input data relative to a particular phase of the 
program and to a particular radial location. As previously 
mentioned, the code handles four input radial locations, and a 
complete set of input data is required at each radial 
location . 

The actual parameters in each file are described in the 
Comments section of the code listing which will follow. They 
are redescribed here to provide the user with a better idea of 
how to specify and use each variable. 

Files 1 through 4 contain airfoil coordinates and other 
information related to the airfoil geometry at each radial 
location. These files are used in computing the flowfield 
about the airfoil. File 1 contains information for the most 
inboard radial station input, with files 2 through 4 
containing input data for progressively more outboard 
stations. 

Card 1 of files 1 through 4 contains the title for the 
Theodorsen transformation input. This title is printed at the 
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beginning of the airfoil section geometry output for the 
particular radial location input. It may typically read 
similarly to, "CLARK-Y COORDINATES AT 30* STATION". 

Card 2 contains five variables: LPT, ITAU, IZ1N, TAU, 

and XZ1N, input in the format (3I5,2F10.0> . 

LPT represents the number of harmonics which the user 
desires in the exponential transformation of the airfoil to 
the circle plane. A value of 50 to 60 is typically used for 
LPT to attain a satisfactory degree of accuracy. 

ITAU is a flag which indicates whether or not the user 
wishes to input the trailing edge angle of the airfoil 
section. A value of 0 indicates that the trailing edge angle 
is not known or will not be input, and any value greater than 
0 indicates that the user will input the section trailing edge 
angle . 

IZ1N is a flag which indicates whether or not the user 
will input the nose singularity coordinates. An IZ1N value of 
0 means that the nose singularity will not be input, and IZ1N 
should be set to any value greater than 0 to input the nose 
singularity coordinates. 

TAU is the trailing edge angle in degrees, to be input 
if ITAU is greater than 0. 

XZ1N and YZ1N are the x- and y-coordinates of the nose 
singularity, to be input if IZ1N is greater than 0. 

Card 3 contains the variables NU and NL, which are the 
number of upper and lower surface coordinates to be input, 
respectively. This card has the format 215. 


89 



' Cards 4 through the last contain the airfoil section 
coordinates, nondimenaionalized with respect to chord. Eight 
values are listed on each line in format 8F10.0. Upper 
surface x-coordinates are all listed first, followed by upper 
surface y-coordinatea, lower surface x-coordinatea, and 
finally lower surface y-coordinatea . / 

■This then completes the input required to map the 
airfoil sections to the circle plane. 

File 8 contains input parameters which are used to 
compute droplet impingement values for the airfoil. The same 
File 8 input data is used for each input radial location. 

Card 1 of file 8 contains the title for that particular 
radial location, and would typically read similarly to "30* 
STATION TRAJECTORY INPUT." This title is printed at the 
beginning of the droplet trajectory section at each radial 
location . 

Card 2 contains the variables IPRINT, MSPL, XFIRST, 

XLAST , and EPS. 

IPRINT is a print option, with a value of 0 used for the 
shortened version of the trajectory summary output, and a 
value of. 1 producing the complete printed output. The 
shortened version will print only the general input to the 
trajectory section, the tangent trajectories, and the 
impingement efficiencies. The trajectory tabulation and 
summary and surface impingement values are not printed when 
IPRINT equals 0. 
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MSPL is a variable which ia uaed to aelect the method 
deaired to calculate the impingement efficiency curve. 

Setting MSPL equal to 0 requeata the program to aelect the 
beat of three methoda poaaible for calculating the curve. An 
MSPL value of 1 cauaea the program to uae a aimple cubic 
apline routine to calculate the curve. An MSPL value of 2 
inatructa the program to uae a cubic apline but to piece in a 
linear fit where the cubic apline faila. An MSPL value of 3 
cauaea a quadratic apline to be uaed, in which the 
coefficienta are aelected to optimize the fit of the curve. 
Setting MSPL equal to 0 haa alwaya produced aatiafactory 
reaulta in all caaea run for thia inveatigation . 

XFIRST ia the x-coordinate after which each atep ia 
checked for droplet impingement on the airfoil. A value of - 
0.05 ia typically uaed for XFIRST. 

XLAST ia the x-coordinate at which droplet trajectory 
calculationa are terminated. Thia value typically rangea from 
0.5 to 1.0, and if the uaer lacks a good feel for the 
impingement characteriatica to be expected for a given caae, a 
value of 1.0 should be used for XLAST to ensure that no 
trajectories are terminated prematurely. 

EPS is the maximum error allowed in the step integration 

- g 

process. A value of l.xlO is generally used for EPS. 

Card 3 contains the variables X0, UD0, Y0, VD0, and FR, 
in format 5F10.5. 

X0 is the initial droplet x-coordinate. Thia ia the x- 
value of the point at which the trajectory inveatigation is 
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begun, and it ia typically given a value of -5.0, representing 
a starting point five chord lengths in front of the airfoil. 

UDO is the initial droplet x-velocity with respect to 
freestream, and is commonly set equal to 1. 

Y0 is the initial droplet y-coordinate . This is the 
starting point for the first trajectory calculation only, and 
its value will depend on the section angle of attack. For a 
positive angle of attack, it is suggested that to maximize the 
chance of initial droplet impingement a slightly negative 
value of Y0 be used, perhaps approximately -0.05. 

VDO is the initial droplet i y- velocity with respect to 
freestream, and is generally equal to 0 in the undisturbed 
flowfield five or more chordlengths in front of the airfoil. 

FR represents Froude number. This parameter may be set 
equal to 0 to ignore the force of gravity on the droplet 
trajectories, and this is commonly done. 

Card 4 contains four variables, DYO, CON, ERRY, and HGT, 
in the format 4F10.6. 

DYO is the increment in Y0 used by the program to vary 
the beginning droplet y-coordinate. Values of DYO are 
commonly on the order of .001 to .01. The code has been set 
up to assign a value for DYO if the user inputs a value of 0 
for this parameter. The internally assigned value will change 
with section angle of attack. If the user inputs a DYO value 
however, this value will be used regardless of the section 
angle of attack. When internally assigned, DYO is given a 
value of 0.003 for section angles of attack less than or equal 
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to 4.6 degrees and is set equal to 0.005 for section angles of 
attack greater than 4.6 degrees. 

CON is a multiplication constant for the impingement 
efficiency curve output. It is generally set equal to 1. 

ERRY is the maximum error acceptable in total collection 
efficiency. A value of 2 <in percent) has commonly been used. 

HGT is a reference height in the y-direction used to 
compute total collection efficiency. If this parameter is set 
equal to 0 the code will calculate and use for HGT the 
projected height of the airfoil section. 

This then completes the droplet trajectory section 

input. 

Files 12 through 15 contain binary flowfield information 
which is output by the flowfield computation section of the 
code when airfoil coordinates (files 1 through 4) are input. 
Since this flowfield information is independent of angle of 
attack and velocity, it is necessary to run the code with 
airfoil coordinates input only once for a particular geometry. 
In subsequent runs using the same geometry, the flowfield 
information calculated and stored in the initial run is read 
and used by the code, thereby eliminating the need for 
recalculation of propeller flowfield information in each run. 
The user must only set up the computer job control language 
(JCL) appropriately to exercise this option. For example, in 
the initial run of a particular geometry, units 12 through 15 
must be set up to have the binary flowfield data written to 
them. In subsequent runs these files must be established as 
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existing data files from which data is to be read. When 
previously computed airfoil flowfield information is used, it 
is not necessary to input airfoil coordinates, and the 
flowfield calculation section of the code will be bypassed. 

File 16 contains input data for the propeller 
performance calculations. This file contains data for the 
entire propeller, so no duplicate files are required for 
different radial locations. Data for a maximum of 15 radial 
locations may be input. 

Card 1 of this file contains the flight and atmospheric 
variables VO, TO, DD, W, TAUM, and ALTUDE in format 6F10.4. 

VO is the freestream velocity in miles per hour. 

TO is the freestream temperature in degrees Fahrenheit. 

DD represents the volume mean droplet diameter of the 
icing cloud in microns. 

W is cloud liquid water content in grams per cubic 

meter . 

TAUM is the icing exposure time in minutes. 

Finally, ALTUDE is the pressure altitude in feet at 
which the icing encounter occurred. 

Card 2 contains the variables NUMBLD, NUMSTA, INCOMP and 
NUMADV in format 415. 

NUMBLD represents the number of propeller blades. 

NUMSTA is the number of radial locations to be input by 
the user (maximum is 15.) 

INCOMP is a compressibility check flag which provides 
the user with the option of obtaining either a compressible or 
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incompressible solution. If set equal to 0, a compressible 
solution is obtained and if set equal to 1, an incompressible 
solution is computed. Recall however that the flowfield 
around the various propeller sections was computed using a 
method which assumes imcompressible flow. 

NUMADV is the number of advance ratios to be input. Set 
equal to 0 if advance ratio is to be calculated using input 
RPM and freestream velocity. NUMADV may have a maximum value 
of 10. 

Card 3 contains the variables RADHUB, BLSET, RPM, 

CPDESI , DESIHP , and DIVIS in format 6F10.5. 

RADHUB is the propeller hub radius in feet. 

BLDSET represents the propeller blade setting in 
degrees. If either CPDESI or DESIHP are assigned values other 
than O, the program will derive a BLDSET value based upon 
these input parameters. 

RPM is the number of propeller revolutions per minute. 

If specified, RPM will be held constant. Otherwise RPM will 
be calculated based on the input freestream velocity and 
advance ratio. Either RPM or advance ratio must be specified 
in the input. 

CPDESI is the design power coefficient. If set equal to 
0, power coefficient will be calculated based upon the input 
value of BLDSET. 

DESIHP is the design horsepower, and if it is input as 
0, a value for DESIHP will be calculated based on the input 
BLDSET value. 
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DIVIS is a convergence aid variable for the design mode 
in which a BLDSET value must be computed. It is generally 
assigned a value between .001 and 1. When BLDSET is 
specified, the DIVIS variable is not employed. 

Card 4 contains the variables CORTIP, RADIUS, STUBLT, 
CDINT, and CORHUB in format 5F10.5. 

CORTIP is the propeller tip chord in feet. 

RADIUS is the propeller radius in feet. 

STUBLT is the propeller stub length in feet. 

CDINT is a blade shank correction factor which takes 
into account interference from the blade shank. If a value 
for CDINT is unknown, a value of 0.567 should be used. 

CORHUB represents the propeller hub chord in feet. 

Card 5 is a title card which denotes the variables on 
the card to follow exactly as shown in the code input 
comments : 

R/RAD= BLDANG= CH0RD= T/C= ALPHA0= CLD= 

VDIS= RC= 

Cards 6 through NUMSTA contain values for these 
variables, in format 8F10.6. 

XPR is the nondimensional radial location to be input. 

BLDANG is the propeller blade angle in degrees of this 
section. Generally this is just the twist distribution built 
into the propeller. 

CHORD is the propeller section chord in feet. 

TC is the section thickness to chord ratio. 
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ALO is the angle of attack for zero lift, in degrees, 
referenced to the longest chordline. 

CLD is the section design lift coefficient. 

VDIS represents the velocity distribution seen by the 
section, taking into account nacelle effects. This parameter 
equals the local velocity divided by that of freestream. 

Finally RC is the section leading edge radius of 
curvature, nondimenaional with respect to section chord. 

The final card(s), in format F10.6, with one value per 
card, are ADVRAT values. 

ADVRAT is the propeller advance ratio, and as many as 10 
values of ADVRAT may be input by the user, with performance 
calculations performed by the code for each advance ratio. If 
NUMADV is set equal to 0, such that advance ratio is to be 
calculated by the code given RPM and freestream velocity, no 
ADVRAT values should be input. 

This completes the propeller performance section input. 

File 17 contains input for options available to the 
user, including choice of aerodynamic coefficient correlation 
to be used as well as print options. 

Card 1 is the title card for the run. This title is 
printed at the beginning of each section of calculations 
within the program and so should be general in nature, such 
as, "ENCOUNTER 2 FROM NEEL AND BRIGHT REPORT." 

Card 2 contains the variables NCD, RADEXT, IREAD, 

IPRCOR , and NSECT in format 15 , F10 . 5 , 315 . 
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NCD is the aerodynamic coefficient correlation selector. 
A value of 1 for NCD calls for the use of Bragg's rime ice AC 

a 

correlation. A value of 2 calls for the use of Gray's glaze 

ice AC. correlation. A value of 3 calls for the use of 
a 

Flemming's AC^ and AC^ correlations for rime and glaze ice. 

RADEXT is the fractional radial extent of icing. 

IREAD is a flag which determines whether or not the code 
will calculate the flowfield at the input radial locations. 

The option is to use previously calculated and stored 
flowfield data. For each radial location input, the code will 
write the flowfield information to a file if it has not been 
computed in a previous run. Set IREAD equal to 0 to calculate 
the flowfield, or set it equal to 1 to use the previously 
computed data. The user must also be sure that the 
corresponding JCL for each of these files is set up correctly, 
i.e., if IREAD equals 0 the JCL must be write-oriented, and if 
IREAD equals 1 the JCL for these files must be read-oriented. 

IPRCOR is a flag which indicates whether or not the user 
wishes to have the airfoil coordinates at each of the four 
input radial locations printed. An IPRCOR value of 0 
instructs the code to not print the coordinates, and an IPRCOR 
value of 1 causes the coordinates to be printed. 

NSECT is a variable which indicates the number of radial 
locations to be input, with a maximum of four possible. 

Cards 3 through 7 contains values for the variable RSTA 
in format FIO.5. 
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RST A is the nondimenalonal radial location at which 


flowfield calculations are to be performed. A maximum of four 
radial locations may be input. These radial locations must 
for the current version of the code be the same as any four of 
the radial locations input in the propeller performance 
section . 

It may be shown that the curves of total and maximum 
local collection efficiencies are relatively smooth functions 
of radial location, so that four points should provide a 
reasonable representation of the variance of these impingement 
efficiencies with radial location. However, it is obvious 
that the fewer points one uses to fit these curves, the less 
accurate they will tend to be. It is recommended therefore 
that whenever possible, no less than the maximum of four 
stations should be input to maximize the available accuracy in 
the impingement efficiency calculations, which then affect the 
propeller performance computations for the iced configuration. 

This then completes the required input to the code. 
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APPENDIX 3 


SAMPLE CASE 

The following output was produced for Encounter 2 at an 
advance ratio of 1.18. The iced calculations were performed 
using Bragg's older, modified correlation. This same format 
is output regardless of the user's choice of correlation. Use 
of the Bragg or Gray correlations will generate output similar 
in appearance to that shown here. When the Flemming 
correlations are used, no flowfield or droplet trajectory 
information is printed out. 

The code requires approximately 570K of computer memory. 
Typical run times are roughly five minutes (CPU time) on an 
IBM 370 and 33 seconds on a CRAY- IS for one advance ratio 
using either the Bragg or Gray correlations. Use of the 
Flemming correlations significantly reduces run time because 
the flowfield and droplet trajectory calculation sections of 
the. code are bypassed, resulting in typical run times for the 
Flemming runs of 8.4 seconds on the IBM 370 and 0.57 seconds 
on the CRAY-1S. 
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Note: Format of files 2-4 identical to format of file 1. Files 2-4 will contain 
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NEEL & BRIGHT DATA, BRAGG CORRELATION 



PROPELLER ICING ANALYSIS PROGRAM 
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STRIP ANALYSIS METHOD 
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DENSITY AT ALTITUDE = 0.0014352 (LBXSECXX2/FTXX4) 

REQUIRED BLADE SETTING = 13.050 (DEGREES) 

REQUIRED HORSEPOWER = 0.000 

REQUIRED POWER C0EFF. = 0.000000 



PROPELLER PERFORMANCE CHARACTERISTICS 
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